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Abstract 

Directionally convex {dcx) ordering is a tool for comparison of dependence 
structure of random vectors that also takes into account the variability of the 
marginal distributions. When extended to random fields it concerns comparison 
of all finite dimensional distributions. Viewing locally finite measures as 
non-negative fields of measure-values indexed by the bounded Borel subsets 
of the space, in this paper we formulate and study the dcx ordering of 
random measures on locally compact spaces. We show that the dcx order 
is preserved under some of the natural operations considered on random 
measures and point processes, such as deterministic displacement of points, 
independent superposition and thinning as well as independent, identically 
distributed marking. Further operations such as position dependent marking 
and displacement of points are shown to preserve the order on Cox point 
processes. We also examine the impact of dcx order on the second moment 
properties, in particular on clustering and on Palm distributions. Comparisons 
of Ripley's functions, pair correlation functions as well as examples seem to 
indicate that point processes higher in dcx order cluster more. 
As the main result, we show that non-negative integral shot-noise fields 
with respect to dcx ordered random measures inherit this ordering from the 
measures. Numerous applications of this result are shown, in particular to 
comparison of various Cox processes and some performance measures of wireless 
networks, in both of which shot-noise fields appear as key ingredients. We also 
mention a few pertinent open questions. 
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1. Introduction 



Point processes (p.p.) have been at the centre of various studies in stochastic 
geometry, both theoretical and applied. Most of the work involving quantitative 
analysis of p.p. have dealt with Poisson p.p.. One of the main reasons being that 
characteristics of Poisson p.p. are amenable to computations and yield nice closed 
form expressions in many cases. Computations have been difficult in great many cases, 
even for Cox (doubly stochastic Poisson) p.p.. 

Comparison of point processes To improve upon this situation, qualitative, com- 
parative studies of p.p. have emerged as useful tools. The first method of comparison 
of p.p. has been coupling or stochastic domination (see [18,20,32]). In our terminology, 
these are known as strong ordering of p.p.. When two p.p. can be coupled, one turns 
out to be a subset of the other. This ordering is very useful for obtaining various bounds 
and proving limit theorems. However, using it one cannot compare two different p.p. 
with same mean measures. An obvious example is an homogeneous Poisson p.p. and 
a stationary Cox p.p. with the same intensity. The question arises of what ordering 
is suitable for such p.p.? This is an important question since it is expected that by 
comparing p.p. of the same intensity one should achieve a tighter bound than by 
coupling. For some more details on strong ordering of p.p. and need for other orders, 
see remarks in [29, Section 5.4 and Section 7.4.2]. 

From convex to dcx order Two random variables X and Y with the same mean 
E(X) = E(Y) can be compared by how "spread out" their distributions are. This 
statistical variability (in statistical ensemble) is captured to a limited extent by the 
variance, but more fully by convex ordering, under which X is less than Y if and 
only if for all convex /, E(/(A)) < E(/(£>)). In multi-dimensions, besides different 
statistical variability of marginal distributions, two random vectors can exhibit dif- 
ferent dependence properties of their coordinates. The most evident example here is 
comparison of the vector composed of several copies of one random variable to a vector 
composed of independent copies sampled from the same distribution. A useful tool 
for comparison of the dependence structure of random vectors with fixed marginals is 
the supermodular order. The dcx order is another integral order (generated by a class 
of dcx functions in the same manner as convex functions generate the convex order) 
that can be seen as a generalization of the supermodular one, which in addition takes 
into account the variability of the marginals (cf [29, Section 3.12]). It can be naturally 
extended to random fields by comparison of all finite dimensional distributions. 
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The dcx order of random measures In this paper we make an obvious further 
extension that consists in dcx ordering of locally finite measures (to which belong 
p.p.) viewed as non-negative fields of measure- values on all bounded subsets of the 
space. We show that the dcx order is preserved under some of the natural operations 
considered on random measures and point processes, such as independent superposition 
and thinning. Also, we examine the impact of dcx order on the second moment 
properties, in particular on clustering, and Palm distributions. 

Integral shot-noise fields Many interesting characteristics of random measures, 
both in the theory and in applications have the form of integrals of some non-negative 
kernels. We call them integral shot-noise fields. For example, many classes of Cox 
p.p., with the most general being Levy based Cox p.p. (cf. [14]), have stochastic 
intensity fields, which are shot-noise fields. They are also key ingredients of the recently 
proposed, so-called "physical" models for wireless networks, as we will explain in what 
follows (see also [1,8,11]). It is thus particularly appealing to study the shot-noise 
fields generated by dcx ordered random measures. 

Since integrals are linear operators on the space of measures, and knowing that a 
linear function of a vector is trivially dcx, it is naturally to expect that the integral 
shot-noise fields with respect to dcx ordered random measures will inherit this ordering 
from the measures. However, this property cannot be concluded immediately from the 
finite dimensional dcx ordering of measures. The formal proof of this fact that is 
the main result of this paper involves some arguments from the theory of integration 
combined with the closure property of dcx order under joint weak convergence and 
convergence in mean. 

Ordering in queueing theory and wireless communications The theory of 
stochastic ordering provides elegant and efficient tools for comparison of random objects 
and is now being used in many fields. In particular in queueing theory context, in [33], 
Ross made a conjecture that replacing a stationary Poisson arrival process in a single 
server queue by a stationary Cox p.p. with the same intensity should increase the 
average customer delay. There have been many variations of these conjectures which 
are now known as Ross-type conjectures. They triggered the interest in comparison 
of queues with similar inputs ( [6,25,31]). The notion of a dcx function was partially 
developed and used in conjunction with the proving of Ross-type conjectures ( [21,22, 
34]). Much earlier to these works, a comparative study of queues motivated by neuron- 
firing models can be found in [16]. Also comparison of variances of point processes and 
fibre processes was studied in [36] and hence it can be considered as a forerunner to our 
article. The applicability of these results has generated sufficient interest in the theory 
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of stochastic ordering as can be seen from the diverse results in the book of Miiller 
and Stoyan ( [29]). As most works on ordering of p.p. were motivated by applications 
to queueing theory, results were primarily focused on one-dimensional point processes. 
An attempt to rectify the lack of work in higher dimensions was made in [24] , where 
comparison results for shot-noise fields of spatial stationary Cox p.p. were given. The 
results of [24] are the starting point of our investigation. 

Our interest in ordering of point processes, and in particular in the shot-noise 
fields they generate, has roots in the analysis of wireless communications, where these 
objects are primarily used to model the so called interference that is the total power 
received from many emitters scattered in the plane or space and sharing the common 
Hertzian medium. According to a new emerging methodology, the interference-aware 
stochastic geometry modeling of wireless communications provides a way of defining 
and computing macroscopic properties of large wireless networks by some averaging 
over all potential random patterns for node locations in an infinite plane and radio 
channel characteristics, in the same way as queuing theory provides averaged response 
times or congestion over all potential arrival patterns within a given parametric class. 
These macroscopic properties will allow one to characterize the key dependencies of 
the network performance characteristics in function of a relatively small number of 
parameters. 

In the above context, Poisson distribution of emitters/receiver/users is often too 
simplistic. Statistics show that the real patterns of users exhibits more clustering 
effects ("hots spots") than observed in an homogeneous Poisson point processes. On 
the other hand, good packet-collision-avoidance mechanisms scheme should create some 
"repulsion" in the pattern of nodes allowed to access simultaneously to the channel. 
This rises questions about the analysis of non-Poisson models, which could be to some 
extent tackled on the ground of the theory of stochastic ordering. Interestingly, we 
shall show that there are certain performance characteristics in wireless networks that 
improve with more variability in the input process. 

The remaining part of the article is organized as follows. In the next section, 
we will present the main definitions and state the main results concerning dcx ordering 
of the integral shot-noise fields. Section [3] will explore the various consequences of 
ordering of random measures. The proofs of the main results are given in Section 0J 
Examples illustrating the use and application of the theorems shall be presented in 
Section \5[ Section [6] will sketch some of the possible applications of results in the 
context of wireless communications. Finally, we conclude with some remarks and 
questions in Section [7] There is an Appendix (Section [8]) containing some properties 
of stochastic orders and their extensions that are used in the paper. 
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2. Definitions and the Main Result 



The order < on R™ shall denote the component-wise partial order, i.e., (x\, . . . , x n ) < 
(yi, ... , y n ) if x t < yi for every i. 

Definition 2.1. • We say that a function / : R d — > R is directionally con- 
vex (dcx) if for every x, y,p,q G K d such that p < x,y < q and x + y = p + g, 

/(*) + /(!/) < f( P ) + /(?). 

• Function / is said to be directionally concave (dcv) if the inequality in the above 
equation is reversed. 

• Function / is said directionally linear (dl) if it is dcx and dcv. 

Function / = . . . , /„) : R d — > M™ is said to be dcx(dcv) if each of its component 
fi is dcx(dcv). Also, we shall abbreviate increasing and dcx by idea; and decreasing 
and dcx by defcx. Similar abbreviations shall be used for dcv functions. Moreover, we 
abbreviate non-negative and idcx by idcx + . 

In the following, let $ denote some class of functions from K d to R. The dimension 
d is assumed to be clear from the context. Unless mentioned, when we state E(/(X)) 
for / e 5 and X a random vector, we assume that the expectation exists, i.e., for each 
random vector X wc consider the sub-class of # for which the expectations exist with 
respect to (w.r.t) X. 

Definition 2.2. • Suppose X and Y are real- valued random vectors of the same 
dimension. Then X is said to be less than Y infi order if E(f(X)) < E(/(Y)) for 
all / G # (for which both expectations are finite). We shall denote it as X <g Y. 

• Suppose {X(s)} s6 5 and {Y(s)} se s arc real-valued random fields, where S is 
an arbitrary index set. We say that {X(s)} <g {Y(s)} if for every n > 1 and 
s u . . . , s n G S, (X( Sl ), . . . , X(s n )) < s (F( Sl ), . . . , Y(s n )). 

In the remaining part of the paper, we will mainly consider ^ to be the class of dcx, 
idcx and idev functions; the negation of these functions give rise to dcv, ddev and ddcx 
orders respectively. If 5 is the class of increasing functions, we shall replace J by st 
(strong) in the above definitions. These are standard notations used in literature. 

As concerns random measures, we shall work in the set-up of [17]. Let E be a 
locally compact, second countable Hausdorff (LCSC) space. Such spaces are polish, 
i.e., complete and separable metric space. Let B(E) be the Borel c-algebra and Bb(E) 
be the a- ring of bounded, Borel subsets (bBs). Let M = M(E) be the space of non- 
negative Radon measures on E. The Borel cr-algebra M is generated by the mappings 
H 1 > m(-B) for all B bBs. A random measure A is a mapping from a probability space 
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(fi, F, P) to (M, M). We shall call a random measure $ a p.p. if $ G N, the subset of 
counting measures in M. Further, we shall say a p.p. <I> is simple if a.s. <f>({x}) < 1 for 
all i eE. Throughout, we shall use A for an arbitrary random measure and $ for a 
p.p.. A random measure A can be viewed as a random field {A(B)}s eBi) (E). With this 
viewpoint and the previously introduced notion of ordering of random fields, we define 
ordering of random measures. 

Definition 2.3. Suppose Ai(-) and A 2 (-) are random measures on E. We say that 
Ai(-) <dcx Aa(-) if for any h, . . . , I n bBs in E, 

(Ai(7i), . . . , Ai(7„)) < dcx (A 2 (/ 1 ),...,A 2 (/„)). (1) 

The definition is similar for other orders, i.e., when $ is the class of idcx/idcv/ddcx/ddcv/st 
functions. 

Definition 2.4. Let S be any set and E a LCSC space. Given a random measure A on 
E and a measurable (in the first variable alone) response function h(x, y) :Ex5-* M + 
where K + denotes the completion of positive real-line with infinity, the (integral) shot- 
noise field is defined as 

V A (y)= f h(x,y)A(dx). (2) 
Je 

With this brief introduction, we are ready to state our key result that will be proved 
in Section |4~T1 

Theorem 2.1. 1. If At <idcx(res P .idcv) A 2j then {VA.i(y)}j/eS <idcx(res P .idcv) 
{V\ 2 (y)} y <=s- 

2. LetE(V Al {y)) < oo, for ally e S,i = 1,2. J/Ai < dcx A 2 , then {V Al {y)} v eS <dcx 

The first part of the above theorem for the one-dimensional marginals of bounded 
shot-noise fields generated by lower semi-continuous response functions is proved in [24] 
for the special case of spatial stationary Cox p.p.. It is conspicuous that we have 
generalized the earlier result to a great extent. This more general result will be used in 
many places in this paper, in particular to prove ordering of independently, identically 
marked p.p. (Proposition 13. 2|) . Ripley's functions (Proposition 13. 4| . Palm measures 
(Proposition I3.5|) . independently marked Cox processes (Proposition I3.7[) . extremal 
shot-noise fields (Proposition 14. ip . Apart form these results, Sections [5] and [6] shall 
amply demonstrate examples and applications that shall need Theorem 12. II 
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3. Ordering of Random Measures and Point Processes 

We shall now give a sufficient condition for random measures to be ordered, namely 
that the condition (JTJ) in Definition 1 2 . 31 needs to be verified only for disjoint bBs. The 
necessity is trivial. This is a much easier condition and will be used many times in the 
remaining part of the paper. 

Proposition 3.1. Suppose Ai(-) and A2(-) are two random measures on E. Then 
Ai(-) <dcx A.2(-) if and only if condition {!]) holds for all mutually disjoint bBs. The 
same results holds true for idcx and idcv order. 

Proof. We need to prove the 'if part alone. We shall prove for dcx order and 
the same argument is valid for / being idcx or idcv. Let condition {T]) be satisfied 
for all mutually disjoint bBs. Let / : R" — » R be dcx function and Bi,... ,B n be 
bBs. We can choose mutually disjoint bBs A\,. . . ,A m such that Bi = \Jj e j i Aj for 
all i. Hence A(B l ) = T,j &Ji H A j)- Now define g : R™ -> as g(xi,...,x m ) = 
QZjgjj x oi ■ • • j Sjej„ x j)- Then g is idl and so / o g is dcx. Moreover, f(A(B\), . . . , 
A(B n )) = f o g(A(A\), . . . , A(vl m )) and thus the result for dcx follows. □ 

3.1. Simple Operations Preserving Order 

Point processes are special cases of random measures and as such will be subject 
to the considered ordering. It is known that each p.p. $ on a LCSC space E can be 
represented as a countable sum $ = EXi of Dirac measures (e x (A) = 1 if x € A and 
otherwise) in such a way that Xi are random elements in E. We shall now show that 
all the three orders dcx, idcx, idcv preserve some simple operations on random measures 
and p.p., as deterministic mapping, independent identically distributed (i.i.d.) thinning 
and independent superposition. 

Let (j> : E — > E' be a measurable mapping to some LCSC space E'. By the image of 
a (random) measure A by <f> we understand A'(-) = A(0~ 1 (-)). Note that the image of 
a p.p. $ by <j) consists in deterministic displacement of all its points by <f). 

Let $ = ^2i£ Xi - By i.i.d. marking of with marks in some LCSC space E', 
we understand a p.p. on the product space E x E', with the usual product Borel 
cr-algebra, defined by $ = J2i £ (xi,Zi), where {Z{\ are i.i.d. random variables (r.v.), 
so called marks, on E'. By i.i.d. thinning of $, we understand $ = J^i ^iS Xi , where 
Zi are i.i.d. 0-1 Bernoulli random variables r.v.. The probability P{Z = 1} is called 
the retention probability. Superposition of p.p. is understood as addition of (counting) 
measures. Measures on Cartesian products of LCSC spaces are always considered with 
their corresponding product Borel a-algebras. 

Proposition 3.2. Suppose Aj,i = 1,2 are random measures and = 1,2 are p.p.. 
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Assume that A\ ^dcx (resp. idcx:idcv) A-2 and $1 "^dcx (resp. idcx: idcv) ^2 ■ 

1. Let A£ fee the image of Ai, i = 1,2, fey some mapping <p : E — > E'. TTien 
Ax <dcx(resp.idcx;idcv) A 2 . ^4s a special case, the same holds true for the dis- 
placement of points of $i 's fey 

2. Let i = 1, 2, be simple p.p. and <I>i, i = 1, 2, fee £/ie corresponding i.i.d. marked 
p.p. with the same distribution of marks. Then 5>i <dcx(resp.idcx-,idcv) 

5. TTien $^ fee i.i.c!. thinning of <&i, i = 1,2, the same retention probability. 

Then Q\ — dcx(resp.idcx;idcv) ^2- 

4- Let Ai and A' 2 fee two random measures such that A^ <dcx(resp.idcx-idcv) A 2 . 
Assume that A^'s are independent of Ai's. Then Ai + A^ <dcx(resp.idcx-,idcv) 
A2 + A' 2 , where + is understood as the addition of measures. 

5. Suppose the random measures are on the product space Ext. Then Ai(E x 
') "Sidcx (resp. idcx; idcv) A 2 (E x •), provided the respective projections are Radon 
measures. 

Proof. 1. The result follows immediately from the Definition 12.31 

2. We shall prove $1 <dcx ^2 and the proof for the other orders is similar. Since E 
is a LCSC space, there exists a null-array of partitions {B n j C E} n >i,j>ij i.e., 
{B n ,j}j>i form a finite partition of E for every n and maxj>i{|i3 nj |} — > as 
n — > 00 where | • | denotes the diameter in any fixed metric (see [17, page 11]). 
For every x € E, let j(n, x) be the unique index such that x 6 B n ji n<!C y Let Z = 
{Z n j} n >i j>i be a family of E'-valued i.i.d. random variables with distribution 
F(-). Define marked p.p. $f = J2x k e*i £ (Xk,z n , Hn , Xk) ) for a = 1,2. We shall 
now verify that the sequences $™'s satisfy the assumption of Lemma [8.21 with 
limits 3>i's respectively. 

Firstly let B 1 ,...,B m C E x E' be bBs and 3 : R m -> R be a continuous 
bounded function. Since B^'s are bounded and $i's are simple, given $i, j = 1,2, 
there exists a.s. N(§i) G N such that for n > AT($j), the indices j(n,Xk) ^ 
i(n, JTj) for X k ^ X u X k ,X t e n (B 1 U . . . U £ TO ). Hence for n > N(^), 
E{g(^{Bi), *?(B m ))|$i) = E(5(*i(Bi), . . . , and in consequence 

Efo(4? (Bt), . . . , $?(B m ))\&i) -> E(5(*i(Bi), . . . , $ i ( J B m ))|$ i ) a.s. as n -> cx>. 
Since g is bounded, by dominated convergence theorem we have that E(g(l > ™(i3 1 ), 
...,*?(B m ))) -> E(g(h(B 1 ),...,^ i (B m ))). Thus (^(Bx),...,^^)) ^ 
(#f(Bi), . . . , $f(B m )). Secondly it is easy to check that for B\ = B' X B", we 
have E($f (Si)) = E($ i ( J B'))- F ( s ") = E($,(Bi)) and hence by an appropriate 
approximation E($f (2?i)) = E($ i (Bi)) for any bBs Bi. 
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Finally for any bBs B c IE x E' and any realization Z =~z = {z n ,j}n>i,j>i, define 
Vf(B) := J E l[(x, 2nj-(n,a;)) £ -B]$i(dac). Since z n jf nr ) is a piecewise constant 
function, l[(x,^ n jV n)X \) G B] is a measurable function in x and so Vfs are 
integral shot-noise fields (as per Definition 12. 4jl indexed by bBs of E x E'. Thus 
from Theorem l2.il we have that for any dcx function /, 

E(/(#?(Bi), . . . , ^(B m ))\Z = z) = E(f(V? (B^, Vf{B m ))) 
< E(/(l?(Bi), . . . , Vi(B m ))) = E(/($J(Bi), . . . , = z) 

Now, taking further expectations we get ( l>™ (£1 ),..., $™(i? m )) < dca; ($%(Bi), 
. . . , $2(Sm))' Since the approximation satisfies the assumption of Lemma I5T21 
the proof follows. 

3. We need to prove E(f(^ 1 (A 1 ),...,^ 1 (A n ))) < E(/($ 2 (^), . . . for 
cfca; (resp. idc; icicv) function / and mutually disjoint Ak, k = 1, . . . , n; cf. Propo- 
sition O Note that given $(A k ) = n k , we have ®(A k ) = ^ where ^ 
are i.i.d. copies of the Bernoulli thinning variable. Thus the result follows from 
the first statement of Lemma 18.31 

4. Using the following fact from [29] : X < dcx (rcsp . idcx . idcv) Y implies X+Z < dcx (rcsp . idcx . ldcv) 
Y + Z provided Z is independent of X and Y one can easily show that Ai + 

Ai <dcx (resp. idcx; idev) A 2 + Ai assuming Ai independent of A 2 . The same 
argument shows that A 2 + A[ < dcx (resp. idcx-. idev) A 2 + A 2 . The result follows 
by the transitivity of the order. 

5. This result follows easily from Lemma 18.21 using an increasing approximation of 
E by bBs. 

□ 

3.2. Impact on Higher Order Properties 

We will state now some results involving ordering of moments of random measures 
and draw some conclusions concerning the so called second order properties. These 
latter ones make it possible to characterize the clustering in p.p.. 

By the n th power of random measure A, we understand a random measure A fe on 
the product space E fc given by A k (Ai x ... x A/.) = Y[j=i^(^-j)- Its expectation, 
a k (-) = E(A fc (-)) is called the kth moment measure. The first moment measure a(-) = 
is called the mean measure. 

Proposition 3.3. Consider random measures Ai <i dcx A 2 . Then K\ <i dcx A 2 and 
< a 2 (')- Moreover, if A% < dcx A 2 then a±(-) = a 2 (-). 
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Proof. By the standard arguments, one can approximate any bBs set Cj, i = 1, . . . , n 
in E fc by increasing unions of rectangles. By Lemma 18.21 and using a similar argument 
about composition of a idcx and idl function as in the proof of Proposition ^. H to prove 
the first statement, it is enough to show the respective inequality for idcx function 
/ : R™ — * R taken of the values of the moment measures on n rectangles in E fc . In this 
context, consider g : R m — > R given by 



g(yi, y m ) = /( n II y j) ' 



where Ji, . . . , J n are fc-element subsets of the set {1, . . . , to}. Note for non-negative 
arguments that if / is idcx then g is idcx. 

The second statement follows easily from the first one by the fact that f(x) = x is 
idcx. For the first moment (mean measure) note that both f(x) — x and f{x) = —x 
are dcx. □ 

We shall explore now the relation between dcx ordering and clustering of points in 
a p.p. One of the most popular functions for the analysis of this effect is the Ripley's 
K function K{r) (reduced second moment function); see [35]. Assume that $ is a 
stationary p.p. on R d with finite intensity A = a(B), where B is a bBs of Lebesgue 
measure 1. Then 

^ (r) = A^ E ( £ mBxAr))-!)), 

where B x (r) is the ball centered at x of radius r and \G\ denotes the Lebesgue measure 
of a bBs G; due to stationarity, the definition does not depend on the choice of G. 

Proposition 3.4. Consider two stationary p.p. i = 1,2, with same finite intensity 
and denote by Ki(r) their Ripley's K functions. If 4>i <d C x $2 then Ki(-) < i^2( - )- 



Proof. Denote I t = E(X)x j6 » 4 nG(*<( S x,( r )) - 1)J, i = 1,2. By the equality of 
mean measures (Proposition 13. 3j) . it is enough to prove that I\ <I%. Note that 7j can 
be written as the value of some shot noise evaluated with respect to 3>f, the second 
product of the p.p.. 

h= l[^ J -eG]l[0<|X fe -X i |<r], 

where denotes the indicator function. Thus, the result follows from Proposition 
and Theorem 12. II 



Another useful characteristic is the pair correlation function defined on R 2 as g(x, y) — 
P2 ^ x '^ , where pk is the kth product intensity, equal (outside the diagonals) to the 



pi(x)pi(y) 
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density of the k th moment measure a k with respect to the Lebesgue measure. 

We avoid discussion on questions such as existence etc. The following result follows 
from Proposition [ 



Corollary 3.1. Consider p.p. such that $1 <dcx 3>2- Then their respective pair 
correlation functions satisfy gi(x,y) < g2(x,y) almost everywhere with respect to the 
Lebesgue measure. 

3.3. Impact on Palm Measures 

For the following definitions and results regarding Palm distributions of random 
measures see [17, Section 10]. 

Definition 3.1. For a fixed measurable / such that < E(J E f(x)A(dx)) < oo, the 
/-mixed Palm version of A, denoted by Af 6 M, is defined as having the distribution 

E(Lf(x)A(dx)l[AGM}) 

P A/ G Af = JE i v r , j „ M G M. 

E{f E f(x)A{dx)) 

In case A (say on the Euclidean space E = R d ) has a density {X(x)} xeR d, we define for 
each x G K d the Palm version A, of A by the formula 

»^ E(A(x)l[A G M]) nr tJ 

P A, G M) = V ^, M G M. 

E(A(x)) 

Palm versions A, can be defined for a general random measure via some Radon- 
Nikodym derivatives. However, we shall state our result for A, as defined above as 
well as for mixed Palm versions Ay in order to avoid the arbitrariness related to the 
non-uniqueness of Radon-Nikodym derivatives. 

Proposition 3.5. Suppose A,-,i = 1,2 are random measures. 

1. If Ai <dcx A 2 then (Ai)f <idcx (A2)/ for any non-negative measurable function f 
such that < f E f(x)a(dx) < 00, where a is the (common) mean measure of Ai, 
i = 1,2. 

2. Suppose that Ai has locally finite mean measure and almost surely (a.s.) locally 
Riemann integrable density \i, i — 1,2. If {Xi(x)} <dcx {^2(2^)}, then Ai <d C x 
A 2 and for every x G R d , (Ai) a < idc x (A 2 ) x - 

Proof. 1. Denote 7j = J E f(x)Ai(dx), i = 1,2. By Proposition 13.31 Ai <d ca; A 2 
implies that the mean measures are equal and thus E(ii) = E(/2). It remains to 
prove 

E( ff (A 1 (S 1 ), . . . , Ai(S n ))7i) < E(g(A 2 (B 1 ), . . .,A 2 (B n ))I 2 ) 
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for idcx function g. This follows from Theorem 1 2 . 1 1 and the fact that h(xo,x) = 
xog(x) : R n+1 — > K is idcx, for non-negative argument o;q. 

2. The first part follows immediately from the second statement of Lemma 18741 For 
the second part, use the same argument about h(xo, x) = xog{x) as above. 

□ 

Remark 3.1. Compared to earlier results where dcx ordering led to dcx ordering, one 
might tend to believe that the loss here (as dcx implies idcx only) is more technical. 
However the following illustrates that it is natural to expect so: consider a Poisson p.p. 
$ and its (deterministic) intensity measure a(-) (i.e., its mean measure a(-) — E($(-)). 
Using the complete independence property of the Poisson p.p. and the fact that each 
dcx function is component-wise convex, one can show that for disjoint bBs Ai, . . . , A n 
and any dcx function /, /(a(Ai), . . . , a{A n )) < E(/($(Ai), . . . , $(A„)). Thus a < dcx 
It is easy to see that a/(-) = a(-) (mixed Palm version of a deterministic measure 
is equal to the original measure). Take f(x) — l[x £ A] for some bBs A. Then 
E($/(A)) = E((<f>(A)) 2 )/a(A) = a(A) + 1 since is a Poisson r.v.. Thus a f (A) < 

E(<5/(A)) disproving atf(A) <dcx ®f(A). Another counterexample involving Poisson- 
Poisson cluster p.p. will be given in Remark 15.21 

3.4. Cox Point Processes 

We will consider now Cox p.p. (see e.g. [35, III 5.2]), known also as doubly stochastic 
Poisson p.p., which constitute a rich class often used to model patterns which exhibit 
more clustering than in Poisson p.p.. 

Recall that a Cox (A) p.p. 4>a on E generated by the random intensity measure 
A(-) on E is defined as having the property that $a conditioned on A(-) is a Poisson 
p.p. with intensity A(-). Note that Cox p.p. may be seen as a result of an operation 
transforming some random (intensity) measure into a point (Cox) p.p.. 

One can easily show that this operation preserves our orders. 

Proposition 3.6. Consider two ordered random measures Ai <dcx (reap. idcx; idev) ^2- 

Then $Ai "^idcx (resp. idcx; idev) ^*A2 ■ 

Proof. Taking a dcx (resp. idcx; idev) function <f>, assuming (by Proposition 13. 1|) 
mutually disjoint bBs Ak, k — 1, . . . , n, using the definition of Cox p.p. and the second 
statement of the Lemma 15731 one shows for i = 1, 2 that that the conditional expectation 

E(c/ ) (<i> M (A 1 ),...^ Ai (A n ))\A l ) 

given the intensity measure Ai is a dcx (resp. idcx; idev) function of (Ai(^4i), . . . , Ai(^4„)) 
The result follows thus from the assumption of the measures Aj being dcx ordered. □ 
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We will show now using Theorem 1 2 . 1 1 that dcx, idcx, idev ordering of Cox intensity 
measures is preserved by independent (not necessarily identically distributed) marking 
and thinning, as well as independent displacement of points of the p.p.. 

By independent marking of p.p. $ on E with marks on some LCSC space E', we 
understand a p.p. $ = Yli e (x i ,z i ) such that given $ = J2i £ xn %i are independent 
random elements in E', with distribution P{Zi £ -|$ = J2i £ xi} = ^ Xi (-) given by 
some probability (mark) kernel F x (-) from E to E'. The fact that F x (-) may depend 
on x (in contrast to i.i.d. marking) is sometimes emphasized by calling $ a "position 
dependent" marking. Independent thinning can be seen as the projection on E of the 
subset $(■, {1}) of the independently marked p.p. $ where the marks Zi £ {0, 1} = E', 
are independent Bernoulli thinning variables Zi = Zi(x), whose distributions may 
be dependent on Xi. Similarly, the projection of an independently marked p.p. $ = 
Y^i £ (xi,Zi) on the space of marks E'; i.e., $(E x •) = J^i £ Zi can be seen as independent 
displacement of points of <& to the space E'. Special examples are i.i.d. shifts of points 
in the Euclidean space, when Zi = xi + Yi, where Yi are i.i.d. 

Proposition 3.7. Suppose = 1,2, are two Cox(Ai) p.p.. Assume that their 

intensity measures are ordered Ai <dcx (resp. idcx- idev) A 2 . Let $>i,i = 1,2 be the 
corresponding independently marked p.p. with the same mark kernel F x (-). Then 

^1 — dcx (resp. idcx; idev) ^2* 

From the above Proposition, the following corollary follows immediately by the last 
statement of Proposition 13. 21 

Corollary 3.2. Independent thinning and displacement of points preserves dcx (resp. 
idcx; idev) order of the intensities of Cox p.p.. 

Proof. (Prop. \57i\ Let $i be Cox (Ai) i = 1, 2 respectively. Assume Ai <dcx(idcx,idcv) 
A 2 . It is known that independent marking of Cox(Aj) p.p. is a Cox(Ai) p.p. with 
intensity measure A, on E x E' given by A^(-) = J E J w l[(x, y) £ ■} F x (dy)Ai(dx); 
cf. [35, Sees 4.2 and 5.2]. Let S be the family of bBs in E x E'; for x £ E and 
bBs C C E x E' consider h(x,C) = J w l[(x,y) £ C]F x (dy). Then the integral shot 
noise V Al (C) = J E h(x, C) A,(dx) satisfies V Al (C) = A;(C) for all bBs C. Thus, by 
Theorem 12. II Ai <d xc (resp. idcx-. idxv) ^2 and the result follows from Proposition 13.61 □ 

If A(-) £ M(R d ) a.s has a density {X(x)} xG ]gd with respect to Lebesgue measure then 
the density is referred to as the intensity field of the Cox p.p., which will be called in 
this case Cox (A) p.p. and denoted by 

It is known that Cox p.p. is over-dispersed with respect to the Poisson p.p., i.e., 
Var($i(i?)) < Var(<& 2 (i?)) where $i,$2 are, respectively, Poisson and Cox p.p. with 
the same mean measure. Hence, it is clear that a Cox p.p. can only be greater in dcx 
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order than a Poisson p.p. with the same mean measure. Indeed, in Section [5] we will 
show several examples when this stronger result holds, namely Cox p.p. that are dcx 
ordered (larger) with respect to the corresponding Poisson p.p., as well as Cox p.p. 
dcx ordered with respect to each other. 

3.5. Alternative Definition of dcx Order 

We viewed a random measure as a random field and have defined ordering from this 
viewpoint. Alternatively, one can consider a random measure as an element of the space 
of Radon measures M and define ordering between two M-valued random elements. 
This can be done once we define what is a dcx function on M. The dcx order can be 
defined on more general spaces; [22] extends the notion of dcx ordering to lattice ordered 
Abelian semigroups with some compatibility conditions between the lattice structure 
and the Abelian structure (LOAS + ). The space M can be equipped with the following 
lattice and algebraic structure. Consider the following partial order: for /x, v e M, we 
say n < v if fi(B) < v{B) for all bBs B in E and addition + v)(B) = /j,(B) + v{B). 
Under this definition, the space M forms a LOAS + as required by [22]. Then one can 
define a directionally convex function on M as in Definition ^. II Call it a dcx 1 function. 
This gives rise to dcx 1 order of random measures analogously to the first part of the 
Definition O 

Now we have two reasonable definitions of ordering of random measures. It is easy to 
see that dcx 1 ordering implies dcx ordering. In light of Example 5.1.7 of [29], existence 
of a counterexample to the converse looks plausible, though we failed in our attempts 
to construct one. However, the result of [3] proves that convex ordering of real valued 
stochastic process {X n }neN implies continuous, convex ordering of the corresponding 
elements of the infinite-dimensional Euclidean spaces M N . This suggests that dcx of 
random measures may imply a dcx 1 * order induced by some subclass of dcx 1 functionals 
of random measures, which are regular in some sense. Leaving this general question 
as an open problem, we remark only that the integral shot-noise fields studied in the 
next section can be seen as some particular class of functionals of random measures, 
which are dcx 1 (in fact linear on M) and regular enough for their means to satisfy the 
required inequality provided the random measures are dcx ordered. It is natural thus 
to have them in the suggested dcx 1 * class. 

Recall also that for strong order of p.p. there is the full equivalence between these 
two definitions, and both imply the possibility of a coupling of the ordered p.p. such 
that the smaller one is a.s. a subset of the greater one; cf [32]. 
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4. Ordering of Shot-Noise Fields 

In this section we will prove Theorem 12.11 concerning dcx ordering of integral shot- 
noise fields, which is the main result of this paper. We will also consider the so called 
extremal shot-noise fields. 

4.1. Integral Shot-Noise Fields 

Usually shot-noise fields are defined for p.p. as the following sum (thus sometimes 
called additive shot-noise fields) V$(y) = J2x„e<i M^ni y) where $ = Yl n £x n an d 
h is a non-negative response function. In definition 12.41 we have made a significant 
but natural generalization of this definition. It is pretty clear as to why we call this 
generalization integral shot-noise field. The extension to unbounded response functions 
is not just a mathematical generalization alone. It shall provide us a simple proof of 
ordering for extremal-shot-noise fields for p.p.. 

Now, we shall prove Theorem 12.11 The proof is inspired by [24] . 

Proof. (Theorem \2.1\ We shall prove the second statement first. The necessary 
modifications for the proof of the first statement shall be indicated later on. 

2. We need to show that (V l (yi), . . . ,V 1 (y m )) < dcx (V 2 ( yi ), . . . , V 2 (y m )) for € 
S, 1 < i < m and V^( ) — Vaj (■), j = 1,2. The proof relies on the construction 

of two sequences of random vectors (V k (yi), . . . , V k (y m )), k = 1, 2 , j = 1, 2 

satisfying the assumptions of Lemma HT2l 

Choose an increasing sequence of compact sets K k , k > 1 in E, such that K k / E. 
Since h is measurable in its first argument, we know that there exists a sequence 
of simple functions h k (-,yi),k £ N such that as k — > oo, hk(-,yi) t M'>2/i) f° r 
1 < i < m. They can be written down explicitly as follows: 

hk(;Vi) = lk1[{x e K k : h(x, y t ) — oo}] 

En— 1.. n — 1 ,. , n , , , . 

^-1[{> G K k : -y- < h(x, yi ) < 

71—1 

for 1 < i < m, where -f k = k2 k . Put P kn = {x € K k : ^1 < h(x lVi ) < |-} and 
P koo = {x 6 K k : h(x, yi) — oo} for 1 < i < m and 1 < n < y k . Note that all P kn 
n = 1, . . . , oo are bBs and the sequence of random vectors we are looking for is 

for j = 1,2. By the definition of integral, it is clear that for j = 1,2 as k — > 
°°) {VjHyi), ■■■iV k {Vm)) T (V 3 (yi), • ■ ■ , V^ym)) a.s. and hence in distribution. 
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By monotone convergence theorem, the expectations, which are finite by the 
assumption, also converge. What remains to prove is that for each k E N, the 
vectors are dcx ordered. 

Fix k G N. Now observe that for j = 1,2, i = 1, . . . , to, Vj! (yi) are increasing 
linear functions of the vectors {^j{I kn ) : n = 1, . . . , -f k , oo), j = 1,2. The 
latter are dcx ordered by the assumptions. And since composition of dcx with 
increasing linear functions is dcx, it follows that (V£(yi), ■ ■ ■ ,V£(y m )) <dcx 
(Vi( yi ),...,V k Hy m )). 

1. For vectors (Vj?(yi), . . . , V£(y m )), k = 1,2...., j = 1,2 defined as above, 
f( v k(yi)>---> v k{ym)) T f(V j (yi),...,V j (y m )) a.s. for / idcx (resp. idev) and 
hence 

E{f{V'(y 1 ),...,V>(y m ))) T E(f(V^ yi ),...,V^y m ))), j = 1,2. The proof is 
complete by noting that E f(V£( yi ), . . . , V£(y m )) < E (f(V^( yi ), . . . ,V k 2 (y m )) 
for all k > 1 and / idea; (resp. zdew). 

□ 

4.2. Extremal Shot-Noise Fields 

We recall now the definition of the extremal shot-noise, first introduced in [13]. 

Definition 4.1. Let S be any set and E a LCSC space. Given a p.p. $ on E and 

a measurable (in the first variable alone) response function h(x, y) : E x S — > R, the 
extremal shot-noise field is defined as 

[/■*(»)= sup{MA^,y)}. (3) 

In order to state our result for extremal shot-noise fields, we shall use the lower 
orthant (lo) order. 

Definition 4.2. Let X and Y be random R d vectors. We say X < to Y if P(X <t)> 
P(Y < t) for every t £ R d . 

On the real line, this is the same as strong order (i.e., when J consists of increasing 
functions) but in higher dimensions it is different. Obviously st order implies lo order 
and examples of random vectors which are ordered in lo but not in st are known; 
see ( [29]). Thus, it is clear that the following proposition is a generalization of the 
corresponding one-dimensional result in [24] where the proof method was similar to 
the proof of the ordering of integral shot-noise fields. We shall give a much simpler 
proof using the already proved result. 
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Proposition 4.1. Let $1 < idc v *2- Then {^* 1 (y)} y es <; {t$ 2 (y)} y es- 

Proof. The probability distribution function of the extremal shot-noise can be ex- 
pressed by the Laplace transform of some corresponding (additive) one as follows. Let 
{xi, . . . , x m } C S and (oi, . . . , a m ) € R m . Then 

P(U( yi ) < a it 1 < i < m) = E(TT l[sup{/i(X n , w ) <<Zi}]) 

n 

= E([][]l^(X n ,y,;) <Oi]) 

i n 

= E(]J]Je logl [' l(x "^^ Ql ]) 

z n 

= E (JJ e -E„-iogi[/i(x„^ I )<a ! ]) 

where U{yi) = — logl[/i(X„,j/i) < o»] is an additive shot-noise with response 
function taking values in [0, oo]. The response function is clearly non- negative and 
measurable. The function f(xi, . . . ,x m ) = er^^ Xi is a ddcx function on (—00,00]. 
The result follows by the first statement of Theorem 12.11 □ 

The extremal shot-noise field can be used to define the Boolean model. Given a 
(generic) random closed set (RACS; see [35, Ch. 6]) G, let h((x,G),y) — l[y G x + G]. 

Definition 4.3. By a Boolean model with the p.p. of germs $ and the typical grain 
G we call the random set C(4>, G) = {y : U^(y) > 0} where $ = Yli £ (x i ,G i ) is i-id- 
marking of <I> with the mark distribution equal to this of G. 

We shall call G a fixed grain if there exists a closed set B such that G = B a.s.. We 
shall demonstrate in Section 16.11 as to how one can obtain comparison results for the 
Boolean model using the results of this section. 

5. Examples of dcx Ordered Measures and Point Processes 

In this section, we shall provide some examples of dcx ordered measures and p.p. 
on the Euclidean space E = M. d . The examples are intended to be illustrative and not 
encyclopaedic. The purpose of the examples is to show that there are dcx ordered p.p. 
as well as demonstrate some methods to prove that two p.p. are dcx ordered. Many 
of the examples seem to indicate that p.p. higher in dcx order cluster more, at least 
for Cox p.p.. 



17 



5.1. Ising-Poisson Cluster Point Processes 

Let {A(s)} se j{£i be a stationary random intensity field. Define a new field, which is 
random but constant in space {A m (s) = A(0)} and deterministic constant field {\h{s) = 
E(A(0))}. Cox(A m ) is known as mixed Poisson p.p. and Cox(X h ) is just the well-known 
homogeneous Poisson p.p.. Denote the random intensity measures of the Cox, mixed 
and homogeneous Poisson p.p., by A, A m and respectively (i.e., A(dx) = X(x)dx, 
etc.) It is proved in [24] that A <dcx A m and when {A(s)} is a conditionally increasing 
field, Ah <dcx A. Recall that a random field {X(s)} is a conditionally increasing field 
if for any k and si, . . . , Sfe G M. d the expectation E(f(X(si))\X(sj) = aj V 2 < j < k) is 
increasing in aj for all increasing /. However, no example of a conditionally increasing 
field was given in [24]. Now we construct one. 

Consider the d-dimensional lattice Z d . Let {X(z)} zeZ d be i.i.d. random variables 
taking values in {+1, —1}. Call {X(z)} a (random) configuration of spins. In order to 
obtain a stationary field consider a random shift of the origin of Z d to U with uniform 
distribution on [0, l] d (U independent of {X(z)}). Let the lattice shifted by U be 
denoted by % d . Pick two numbers /i2 < Mi- For s G M. d , define A(s) = /j,il[X(s) = 
1] + H2l[X(s) = — 1] where s represents the unique "lower left" point in 7L d nearest to s. 
The intensity field is clearly stationary. We shall now show that {A(s)} is conditionally 
increasing. Note that 

f(X(s)) = l[x(s) = l](.f( Ml ) - f( P2 )) + /M (4) 

From Theorem 1.2.15 of [29], it is sufficient to show the conditional increasing property 
conditioned on U, the random origin of the lattice 7L d . Hence it is enough for the Ising 
model to possess the following property: 

P(X( Zl ) = l\X(z 2 ) = -l,X( Zj ) = aj , j = 3, . . . , k) 
< P(X( Zl ) = l\X(z 2 ) = l,X( Zj ) = aj , j = 3, . . . , k), 

where G {+1,-1} and Zi G 7L d ,i = 1, . . . , k. This follows easily from the fact that 
the spins are i.i.d. 

We call the Cox p.p. generated by the above conditionally increasing field {A(s)} 
the Ising-Poisson cluster p.p. By the arguments presented in [24], it is dcx larger 
than the homogeneous Poisson p.p. with the same intensity. Note that intuitively the 
Ising-Poisson cluster p.p. "clusters" its points more than a homogeneous Poisson p.p. 
In what follows, we will see more examples of cluster (Cox) p.p. which are dcx larger 
than the corresponding homogeneous Poisson p.p.. 
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5.2. Levy Based Cox Point Processes (LCPs) 

This class of p.p. is being introduced in [14]. One can find many examples of LCPs 
in the above mentioned paper. In simple terms, a LCP is a p.p. whose intensity held 
is an integral shot-noise field of a Levy basis. A random measure L e M(]R d ) is said to 
be a non-negative Levy basis if 

• for any sequence {A n } of disjoint, bBs of M d , L(A n ) are independent random 
variables (complete independence) and L([J A n ) — ^2 L(A n ) a.s. provided UA n 
is also a bBs of R d . 

• for every bBs A of M d , L(A) is infinitely divisible. 

We shall consider only non- negative Levy bases, even though there exist signed Levy 
bases too (see [14]). Hence, we shall omit the reference to non-negativity in future. 
A Cox p.p. <E> is said to be a LCP, if its intensity field is of the form 



where L is a Levy basis and the kernel k is a non-negative function such that k(x,y) 
is a.s. integrable with respect to L and k(.,y) is integrable with respect to Lebesgue 
measure. In [14] the response function k and the Levy basis L is chosen such that 
f B X(y) dy < oo a.s. for all bBs B, for which a sufficient condition is J B E(A(y)) dy < oo. 
In our considerations, in order to be able to use Lemma l8.4l we will require that X(y) 
is a.s. locally Riemann integrable. 

Remark 5.1. Note that a sufficient condition for this is that X(y) is a.s. continuous, 
for which, in turn, it is enough to assume that k is continuous in its second argument 
and that for all x G K d , there exist B x (e x ), e x > such that J Md sup zgB;c( - ex ) k(z, y) a(dx) 
< oo for all y, where a(B) = E(L(B)), the mean measures of the Levy bases; (cf [1]). 

Lemma 5.1. Let L\ and Li be Levy bases with mean measure on. Let = 1,2 be 
LCPs with Levy bases Li,i — 1,2 respectively. 

1- Li <dcx{res P .idcx-idcv) L 2 if and only if L X (A) < C x(resp.icx;icv) L 2 {A) for all bBs 
A ofM. d , where cx,icx,icv stands, respectively for convex, increasing convex and 
increasing concave. 

2- If Li <dcx(res P .idcx-idcv) ^2, then^i <dcx (resp. idcx-,idcv) $2 provided the intensity 
fields Xi(y) of LCP $i is a.s. locally Riemann integrable with these integrals, in 
case of dcx, having finite means. 

3- cti <dcx Li- 
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Proof. The first part is due to Proposition ^ . 1 1 and the complete independence prop- 
erty of Levy bases. As for the second part, it is a simple consequence of Theorem 12. 1[ 
Lemma 18.41 and Proposition 13.61 The third part follows from complete independence 
and Jensen's inequality. □ 

We shall now give some examples of dcx ordered Levy basis. 

Example 5.1. Let {x{\ be a locally finite deterministic configuration of points in E d . 
Let {Xf}i>i,j — 1, 2 be i.i.d sequence of infinite divisible random variables such that 
Xl < cx X\. (For example, X\ can be sum of two independent exponential r.v. with 
mean 1/2 and X^be an exponential r.v. with mean 1.) Define the Levy bases as 
follows: 

LM)=Y, x3 i> 

Xi£A 

where A is a bBs of R d and j — 1,2. By Lemma \5. II and the fact that X\ < cx X\ it 
follows that L\ <dcx L*2- 

Example 5.2. Let $ = X)i be an homogeneous Poisson p.p. on R d inde- 

pendently marked by random variables {Zi} with mean Ao- Consider two random 
measures Ai = J2(x- z )ei> ^o £ *i and A2 = J2( x - z )ei Zi£ Xi - Note that Li, i — 1,2 are 
Levy basis. By Lemma 15.11 and the fact that Ao < CI Zi, conditioning on the number 
of points and using the same argument as in the proof of the second statement of 
Proposition [321 one can prove that Ai <dcx A2. 

5.3. Poisson-Poisson Cluster Point Processes 

By Poisson-Poisson cluster p.p., we understand a LCP with the Levy basis being a 
Poisson p.p. This class deserves a separate mention due to the generality of the ordering 
results that are possible. For rest of the section, assume that h(x) is a non-negative 
measurable function such that J Rd h{x)dx = Xq < 00. 

We shall now give an example of a parametric family of dcx ordered Poisson-Poisson 
cluster p.p.. Fix A > 0. Let $ c , c > be a family of homogeneous Poisson p.p. on Mr 
of intensity cA. Let a non-negative function h : R d x M. d — > K be given and consider 
a family of shot noise fields X c (y) = J R d {h{x,y)/c) & c (dx), which are assumed a.s. 
locally Riemann integrable with J B E(A c (y)) dy < 00 for bBs B. 

Proposition 5.1. The family of shot-noise fields {X c (y)} yGR d is decreasing in dcx, 
i.e., for < c± < c 2 we have {A C2 (y)} < dcx {A Cl (y)}. Consequently Cox(X C2 )< dc x 
Cox(X Cl ). 

Proof. Note that {A c (a:)} can be seen as a shot-noise field generated by the response 
function h and the Levy basis L c — (l/c)Q\ c . By Lemma \5. II and Theorem 12.11 it is 
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enough to prove that L C2 (A) < cx L Cl (A) for A bBs and C2 > c\ > 0. 

Since, X < cx Y implies that aX < cx aY for all scalars a > 0, it suffices to prove that 
L C a(A) < cx L a (A) for A bBs and c > 1, a > 0. This essentially boils down to proving 
that N ca < cx cN a , c > l,a > 0, where N a stands for a Poisson r.v. with mean a. 
Let {^™}i<j< n and {^ n }i<i< n , n > 1 be i.i.d. sequences of Bernoulli r.v's with 
probability of success ca/n and a/n, respectively, with n > ca. Let X n — Y^i=\ an d 
Y n = Y^i=i c ^i™- It is well known that X n , Y n converge weakly to iV ca , N a respectively, 
as n — > oo. As convex order preserves weak convergence, we need to only prove 
that X n < cx cY n . By the independence of summands, it is enough to prove that 
Xf < cx cY™, which we shall do in what follows. Let / be a convex and differentiable 
function. Define g(c) := E f(X") — E J{cY- 1 ) = £{c(/(l) - /(0)) - /(c) +/(0)}. Note 
that g(l) = 0. Hence, our proof is complete if we show that g is decreasing in c > 1. 
Indeed, 

g'(c) = £{(/(!) _/(0))-/'(c)} 

= -{/'(&) -/'(c)}<0, (6<c) 
n 

where 6 G (0, 1) by mean-value theorem and /' is increasing due to convexity. □ 

Poisson-Poisson cluster p.p. can be also dcx compared to a homogeneous Poisson 
p.p.. Let $ and $' be homogeneous Poisson p.p. with intensities A < oo and A x Ao 
respectively. Define fi(y) = J2x e* H^i — y)- Let $" be Cox(/z(x)). 

Proposition 5.2. Let {m(2/)} ^ e as above. Assume that (i(y) is a.s. locally 

Riemann integrable and E(/i(y)) = E(/i(0)) < oo. Then <<jca: 

Proof. By the last statement of Lemma [5.11 we have Ada: <dca; &(dx). Note that 
A x Ao = J Rd h(x — y) Xdx and thus by the second statement of Theorem 12.11 (note 
the assumption E(/i(y)) < oo) {A x A } <d cx {l J '(y)}i where the dcx smaller field 
is a deterministic, constant. The result follows now from the second statement of 
Lemma 18.41 by assumption that fJ.(y) is a.s. Riemann integrable and observing that 
E(J A fx(y) dy) = E(/i(0)) f A dy < oo for all bBs A. □ 

Remark 5.2. Consider Poisson p.p. <&' and Cox(/i) as in Proposition l5.2l It is known 
that the Palm version (given a point at the origin) of can be constructed taking 
+ So- By [27, Proposition 2], analogously, Palm version of Cox(/x) can be taken as 
Cox(/i) + £o + where $" is an independent of Cox(/z) Poisson p.p. with intensity 
h{y ~~ where ^ is sampled from the distribution h(dx)/ J h(y)dy. This shows that 
one cannot expect dcx ordering of the Palm versions of $' and Cox(/i). 
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5.4. Log Cox Point Processes 

This class of p.p. are defined by the logarithm of their intensity fields. 

An extension of LCP studied in [14] is Log-Levy driven Cox process (LLCPs). Under 
the notation of the previous subsection, a p.p. $ is said to be a LLCP if its intensity 
field is given by 



[14] allows for negative kernels and signed Levy measures but they do not fit into our 
framework. Suppose that L± <idcx L 2 , then $1 <idcx ^2 where <E>i,i = 1,2 are the 
respective LLCPs of Li, i = 1,2 with kernel .). These are simple consequences of 
Theorem 12.11 and the exponential function being icx. 

Another class is the Log-Gaussian Cox process (LGCPs)(see [26]). A p.p. $ is 
said to be a LGCP if its intensity field is X(y) = exp{A(y)} where {X(y)} is a 
Gaussian random field. Suppose {Xi(y)},i = 1,2 are two Gaussian random fields, 
then {Xi(y)} < ldcx {X 2 (y)} if and only if E(Xi(y)) < E(A 2 (y)) for all y e K rf and 
cov(Ai(yi), Xi(y 2 )) < cov(A 2 (yi), A 2 (y 2 )) for all 1/1,2/2 G From the composition 
rules of idcx order, it is clear that idcx ordering of Gaussian random fields implies 
idcx ordering of the corresponding LGCPs. An example of parametric dcx ordered 
Gaussian random field is given in [24, Sec 4]. 

5.5. Generalized Shot Noise Cox Processes (GNSCPs) 

This class of Cox p.p. was first introduced and its various statistics were studied 
in [28]. In simple terms, these are Cox p.p. whose random intensity field is a shot- 
noise field of a p.p. We say a Cox p.p. is GNSCP if the random intensity field 
{A(y)}j, gR d driving the Cox p.p. is of the following form : A(y) = . Jjk^ (cj, y) 
where (c,-, bj, 7,-) £ $, a p.p. on R d x (0, 00) x (0, 00). Also we impose the following 
condition on the kernel k : kf )j (cj,y) — kl{ - c ^ h i^ y l h ^ where ki(cj, .) is a density with 
respect to the Lebesgue measure on R d . We shall denote the GNSCP driven by $ as 
$ G . This class includes various known p.p. such as Neyman-Scott p.p., Thomas p.p., 
Matern Cluster p.p. among others. The case when b/s are constants and {(cj,7j)} is 
a Poisson p.p. is called as Shot Noise Cox process (See [27]). Shot Noise Cox process 
are also LCPs. Suppose two p.p. $1 < dcx (icsp. idcx- idcv) $2, then from Theorem 12.1) 
we infer that $f < dcx (resp . idcx . idcv ^ 

5.6. Ginibre-Radii Like Point Process 

Let { < E > i}i>o be an i.i.d. family of p.p. on R + . So, the points of each p.p. <!>i can 
be sequenced based on their distance from the origin. Let $ be the p.p. formed by 
picking the ith point of <&i for i > 1. We shall from now on abbreviate $([0, b]) by $(6) 
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for b > and similarly for other p.p. used. Note the following representation for $(&) 
and $o(&): 

fc>i fc>i 

Let 

m rn 

<F"(&) = ]T l[*fc(6) > k] ; = l[*o(6) > *]• 

fc>i fe>i 

By Lorentz's inequality (see [29, Th. 3.9.8]), it follows that (<& x (&),..., $ m (6)) < sm 
($o(6), . . . , $ o(fr))j where sm stands for supermodular (sec [29, § 3.9]). Define the 
/ : N m — > R as follows : /(n x , . . . , n m ) = ^Zfe>i l[ n fc > k]. It is easy to verify that 
both / and — / are sm and f(n A m) < /(n),/(m) < f(n V m). In consequence 
g o / is sm provided 3 is ca; and E(g($ m (6))) = E(g o /($ x (&), . . . , $ m (6))) < E(g o 
/($o(6),-..,*o(&))) = E(s($^(6))). Hence $ m (6) < C2; $^ l (6) and using Lemma EH 
we get that $(&) < cx ®o{b). To complete the proof $ <dca: $0, one would require 
a multi-variate generalization of Lorentz's inequality which we have been unable to 
prove. 

We shall now explain the reasons for considering the above p.p. <I>. If we assume 
that $i above are Poisson, then 4> is know to be a representation of the p.p. of the 
squared radii |<&g| 2 = {|^«| 2 : X n G $g} of the Ginibre process $g (see [4,19]). It 
has been observed in simulations that this determinental p.p. exhibits less clustering 
than the homogeneous Poisson p.p. Our result can be seen as a first step towards a 
formal statement of this property. 

6. Applications to Wireless Communication Networks 

From the point of view of applications of our main result, what remains is examples 
of interesting dcx functions. In what follows, we will provide such functions arising 
in the context of wireless networks. In many of the models we have assumed ordered 
point processes with i.i.d. marks. However due to Propnosition 13.71 the results hold 
for independently marked Cox p.p. provided the respective intensity measures are 
ordered. 

6.1. Coverage Process with Independent Grains 

The Boolean model C($,G) defined earlier (see Definition 14. 3|) is the main object 
of analysis in the theory of Coverage processes (see [12]). The percolation properties 
of the Boolean model has been studied in [23] while the connectivity properties of the 
Boolean model has been studied in [30] . For $ as in the Definition 14.31 of the Boolean 
model, denote by V(y) = J2(x- G )e* l[f ^ Xi + @i] the number of grains covering 



23 



y £ M. d . Denote by ^(sx, s n ) the joint probability generating functional (p. g.f) of the 
number of grains covering locations y\, . . . , y n G R d ip(si, . . . , s n ) — EM}" =1 sJ^ Vl ^ , 
Sj > 0, j — 1, . . . , n. Note that the function g(v\, . . . , w„) = , Sj 3 is idea; when Sj > 1 
for all j = 1, . . . , n and is ddcx when < s j < 1 for all j . 

Thus the following result follows immediately from Theorem 12. 1[ Proposition 13.21 
and Proposition 13.71 

Corollary 6.1. Let $i, i = 1,2 &e a simple p.p. (of germs) on M. d . Consider the cor- 
responding Boolean models with the typical grain G and, as above, denote the respective 
coverage number fields by {Vi(y)} and and their p. g.f by ipi. J/$i <dcx{reap.idcx;idcv) $2 
then {Vi(y)} <<fc x (resp. idex; idev) {Vsiy)}) with the result for dcx holding provided 
E(Vi(y)) < oo for all y. In particular, if ^ < ldcx $ 2 then ^{V^yf) < E(V 2 (y) p ) 

for all [3 > 1. <idcx {resp. ddcx) *2 then 1pl(si, . . . , S n ) < 1p2{si, ■■■ , S„) for Sj > 1 

f resp. Sj < 1) j = 1, . . . , n. 

Note that 1 — ?/>(0, . . . , 0) represents the expected coverage measure, i.e., the probability 
whether the locations y%, . . . , y n are covered by at least one grain. In [12, Section 3.8] 
it is shown that expected one-point coverage (or volume fraction in case of stationary 
p.p.) for a stationary Cox p.p. and some clustered p.p. is lower than that of a 
stationary, homogeneous Poisson p.p.. 

Coverage processes arise in various applications. In particular, in wireless commu- 
nications the points of the p.p. (germs) usually represent locations of antennas and 
their grains the respective communication regions. In this context V{y) is the number 
of antennas covering the point y and the coverage measure is the indicator that at least 
one of them is able to reach y. The application of the Boolean model to the modeling 
of wireless communications dates back to the article of Gilbert [10] in 1961. 

6.2. Random Geometric Graphs (RGGs) 

This class of graphs has increasingly found applications in spatial networks. For 
a detailed study of these graphs, see [30]. A random geometric graph is defined as 
a graph with $ as the vertex set and the edge-set E = {{Xi,Xj} : \X{ — Xj\ < r}. 
Clearly this is related to the Boolean model defined in the previous subsection. One 
of the objects of interest in a RGG is the typical degree. Under the notation of the 
previous subsection, the typical degree (deg(&, G)) for a RGG formed by a stationary 
p.p. $ and grain distribution G is deg(tf>,G) — xTK\ S y : XjSfl? 

l[Xi £ A]l[Xt ? 

Xj]l[(Xi + d) n (Xj + Gj) i= 0], where A is a bBs. If G = B (r),r > 0, then 
E(deg(&, G)) = K(r) is the Ripley's K function defined in Section I3~2l The following 
result follows easily from Theorem 12. 11 Proposition 13.31 and Proposition 13. 71 

Corollary 6.2. Suppose that simple p.p. $i <dcx &2, then deg(§\,G) <idcx deg($2,G) . 
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6.3. Interference in Wireless Communications 

The Boolean model is not sufficient for analyzing wireless networks as it ignores the 
fact that in radio communications signal received from one particular transmitter is 
jammed by the signals received from the other transmitters. According to information 
theory as well as existing technology, the quality of a given radio communication link is 
determined by the so called signal to interference and noise ratio (SINR) at the receiver 
of this link, a mathematical point of view, the interference in the above considerations 
is just the sum of the powers of the signals received from all transmitters (perhaps 
except own transmitter (s)). It is then the shot-noise field of received powers that plays 
important role in determining the connectivity and the capacity of the network in a 
broad sense. The foundations of the theory of SINR coverage processes are quite recent 
(see [1,2,8, 11]). In what follows, we shall study the impact of structure of the p.p. of 
interferers on given radio links. 

Consider a set of n emitters {x{\ and n receivers {yi}. Suppose that the signal 
received by yi from Xk is Ski- These {Sik} are assumed to be independent. The 
assumption of independence is due to the phenomenon of fading. Let the set of addi- 
tional interferers be modeled by a i.i.d. marked p.p. $ = £(x j .(z 1 ,....z n )j independent 
of {Sik}, where Z™ is the power received by the receiver yi from the interferer located 
at Xj. Denote the background noise random variable by W. 

We say that the signal from Xi is successfully received by yi if Sa/(W + Ii + Vi) > T 
where Ii = J2k^i &ki and Vi = ■ Zj is the interference received at yi from the set 
of other emitters {xk : k ^= i} and interferers in $, respectively, and T > is some 
(assume constant) required SINR threshold. If we denote by p, the probability of 
successful reception of signals from each Xi to yi, then 

p = P(Su> (W + Ii + Vi)T V* = l,...,n) 

= ni[Fn(T(W + h + V % ))), (5) 

i 

where Fa(s) — P(Su > s) and the second equality is due to independence. Given 
{Ii : i = 1, . . . , n} and W, the expression under expectation in ^ can be viewed as a 
function of the value of the shot-noise vector (Pi, ... , V n ) evaluated with respect to 
Theorem 12.11 and Proposition 13 . 71 implies the following result concerning the impact of 
the structure of the set of interferers on p. 

Corollary 6.3. Consider emitters {xi}, receivers {y}i, powers {Ski} o,s above. Let 
u = 1,2 be two simple marked p.p. of interferers. Denote by p u , u — 1, 2 the prob- 
ability of successful reception given by (5p in the model with the set of interferers <!>„ . 
Assume the product of tail distribution functions of the received powers Y\7=i Fa( s i) 
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be dcx. 7/$i <ddcx $2 then pi <P2- 

It is quite natural to assume ddcx Yli=i^ii( s i)- For example the constant emitted 
power P, omni-directional path- loss function l(r) and Rayleigh jading in the radio 
channel implies Ski = PHki/l(\xk — Ui\), where \Hki\ are i.i.d. exponential random 
variables with mean 1. In this case 11™= 1 Fii( s i) is ddcx. (* J 

7. Conclusions and Open Questions 

To the best of our knowledge, this is the first study of dcx ordering of random 
measures and p.p.. We have defined the dcx order and characterized it by finite 
dimensional distributions of the measure values on disjoint bBs of the space. As the 
main result, we have proved that the integrals of some non-negative kernels with respect 
to dcx ordered random measures inherit this ordering from the measures. This was 
shown to be a very useful tool in study of many particular characteristics of random 
measures and in the construction and analysis of stochastic models. 

In this paper, we have also left several open questions. Here we briefly summarize 
them. 

• Our dcx order is defined via finite dimensional distributions of random measures. 
This makes the verification of dcx order more easy but requires additional work 
when studying functionals, which cannot be explicitly expressed in terms of the 
values of the measure on some finite collection of bBs as, e.g., an integral of the 
measure. Considering a dcx 1 * order on the space of measures could facilitate the 
former task. However, the precise regularity conditions of the dcx 1 * functional 
on the space of measures which would guarantee the equivalence between these 
two approaches are not known (cf Section 13. 5p . 

• Comparisons of Ripley's functions (see Proposition 13. 4j) and pair correlation 
functions (Corollary 13. 1|) seem to indicate that the higher in dcx order processes 
cluster their points more. We have shown examples of p.p., which are larger 
than Poisson one, namely Cox p.p., which indeed exhibit more clustering than 
in Poisson p.p.. It would be interesting to show examples of p.p. which are dcx 
smaller than Poisson one, and which exhibit less clustering than it. Matern "hard 



(*) Recently in [9], under the assumption of Rayleigh fading, direct analytical methods have 
been used to compare the probability of successful reception in Poisson p.p. and a class of Poisson- 
Poisson cluster p.p. known as Neyman-Scott p.p. for both stationary and Palm versions. These 
results relay on explicit expressions for this probability known in the considered cases. Further, it is 
shown that for a certain choice of parameters, Palm version of the Poisson-Poisson cluster p.p. has a 
worser probability of successful reception than the Poisson p.p.. In our terminology, it simply means 
that the corresponding Palm versions aren't ddcx ordered as the connectivity probability is a ddcx 
function (Eqn. [5} of the integral shot-noise fields of the corresponding Palm versions. This strengthens 
Remark 1 5. 2 1 by showing that idcx ordering of Palm versions is the best one can obtain in full generality. 
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core" p.p. and Ginibre p.p. are some natural candidates for this. 

• We have studied dcx order that takes into account the dependence structure and 
the variability of the marginals or random measures. It seems plausible to study 
in a similar manner other orders such as convex, component-wise convex order 
etc. Note however that the supermodular order does not seem to be a reasonable 
one in the context of random measures. The reason is that it allows to compare 
only measures with the same finite dimensional distributions, and thus a Poisson 
p.p. can only be (trivially) compared in this order to itself. Indeed, Poisson finite 
dimensional distributions imply total independence property and thus uniquely 
characterize Poisson p.p. (cf [7, Lemma 2.3.1]). 

8. Appendix 

In order to make the paper more self-contained, we shall recall now some basic 
results on stochastic orders used in the main stream of the paper. The following two 
lemmas can be found in [29, Chapter 3]. 

Lemma 8.1. 1. A twice differentiable function f is directionally convex if and only 
d 2 

if o p. f{x) > 0, for all x, 1 < I, j < n. 

2. The stochastic order relation <dcx is generated by infinitely differentiable dcx 
functions. 

Due to the above lemma, at some places we only prove that two random vectors are 
ordered with respect to twice differentiable dcx functions and conclude that they are 
dcx ordered. 

We denote by —> convergence in distribution (weak convergence). 

Lemma 8.2. Let (X^ : k = 1, . . .) and (Y^ : k = 1, . . .) be sequences of random 
vectors. Suppose X^ < dcx for all k e N. If X™ ^ X and Y^ ^ Y and if 

moreover E(X^) -» E(X) and E(Y^) -» E(Y), then X < dcx Y. 

The following result is from Lemmas 2.17 and 2.18 of [21]. 

Lemma 8.3. 1. For i = 1, . . . , m let (Sj : j = 1, . . .) be independent sequences of 
i.i.d. non-negative random variables. Suppose f is dcx (resp. idcx; idev), then 
g(n 1} n m ) = E(/(^"i 1 Sj, Y!j=i S f)) is also dcx ( res P- idcx : idev). 

2. Let Ni,i — 1, . . . , k denote k mutually independent Poisson r.v. where the mean 
of Ni is A,. If <p '■ N fe — > R is dcx {resp. idcx; idev), then g(Xi, . . . , Xk) = 
E((f>(Ni, . . . , iVfc)) is also dcx (resp. idcx; idev). 
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The first part of the following lemma is an easy extension of the one-dimensional version 
in [21]. The second part, which we prove in what follows, is a further extension of it. 

Lemma 8.4. Suppose {X(s)} seR d and {Y(s)} s£ s.d are two non-negative real-valued 
and a.s. locally Riemann integrable random fields. For some n > 1 and disjoint bBs 
Ii,...,I n denote J x = J T X(s)ds, J Y = Jj. Y(s)ds. 

1- If{X{s)} ^idcx(resp.idcv) then (Jx ) • • • > Jx) — idcx (resp . idev) (Jy, ■ ■ ■ , Jy)- 

for any n and for any I\, . . . , I n disjoint bBs. 

2. Suppose further that E(J A X(x)dx) < oo for all bBs A in R rf and similarly for 
{Y(x)}. If {X(x)} < dcx (dcv){Y(x)}, then (J x ,..., J x ) < dcx{dcv) (J Y , Jy). 

Proof. (2) We shall prove for d — 1 and as can be seen from the proof, the 
generalization is fairly straightforward. 

We need to prove that (Jj X(s)ds, ■ ■ ■ , Jj X(s)ds) < dcx (/j Y(s)ds, . . . , f T Y(s)ds), 
for Zj, i = 1, . . . , n disjoint bBs. We shall give an approximation satisfying the assump- 
tions of Lemma 18.21 Let = [a,, &,]; a,, 6j G R, i = 1, ...,n. Let {{t l m j)i<j<k m , i = 
l,...,rt} be the sequences of mth nested partition of each interval. The middle 
Riemann sum can be given as follows : X m (Ii) = Yl,j x ( t% mj)iKn(j+i) _ *toj)>* = 
1, . . . , n, k £ N and similarly for Y(x). These are the variables satisfying the approxi- 
mation as in Lemma 18.21 As X(s) is Riemann integrable, 

(X m (I 1 ),...,X m (I n ))^(J x ,...,J x ) 

a.s. and hence in distribution. It is also clear the middle Riemann sums of X(-) and 
Y(-) are ordered. What remains to prove is that EI m (/i) — > E J\. In the last term, 
by Fubini, we can interchange the expectation and integral and hence it suffices to 
prove EX m (Ii) — * J EX(s)ds. Our assumption implies that this is true. □ 
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